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We introduce a method, based on an exact calculation of the effective action at large N, to bridge 
the gap between mean field theory and renormalization in complex systems. We apply it to a d- 
dimensional manifold in a random potential for large embedding space dimension TV. This yields 
a functional renormalization group equation valid for any d, which contains both the 0(e = 4 — d) 
results of Balents-Fisher and some of the non-trivial results of the Mezard-Parisi solution thus 
shedding light on both. Corrections are computed at order 0(1/N). Applications to the problems 
of KPZ, random field and mode coupling in glasses are mentioned. 



The random manifold problem, i.e. the behavior of an 
elastic interface in a random potential is important for 
many experimental systems and still offers a consider- 
able theoretical challenge 0. It is the simplest exam- 
ple of a class of disordered systems, including random 
field magnets, where the so called dimensional reduction 
H renders conventional perturbation theory trivial and 
useless. It also provides powerful analogies, via mode 
coupling theory, to complex systems such as structural 
glasses [[[) . Two analytical approaches have been devised 
so far, in limits where the problem appears to simplify 
while remaining non-trivial: the functional renormaliza- 
tion group (FRG) H [|, and the (mean field) replica 
gaussian variational method (GVM) together with 
their dynamical versions 0, ||. The FRG is hoped 
to be controlled for small e = 4 — d, where d is the in- 
ternal dimension of the interface (parameterized by a TV 
component vector u(x) in the embedding space). It fol- 
lows the second cumulant of the random potential R(u) 
under coarse graining, which becomes non-analytic at 
T = beyond the Larkin scale. The GVM approxi- 
mates the replica measure by a replica symmetry broken 
(RSB) gaussian, equivalently, the Gibbs measure for u as 
a random superposition of gaussians H , and is claimed 
to be exact for N — oo. Computing next order correc- 
tions is fraught with difficulties |l3| , and it is still unclear 
in which sense both methods describe small but finite 
e or 1/N. The GVM for instance predicts a transition 
for d = which must disappear for any finite N. We 
found recently [l2| [p| that higher loop FRG equations 
for R(u) at u ^ contain non-trivial, potentially ambigu- 
ous "anomalous terms" involving the non-analytic struc- 
ture of R(u) at u = 0. Although a solution was found to 
two loops [|l2| [ll| , the many loop structure remains mys- 
terious. While both methods circumvent dimensional re- 
duction by providing a non-perturbative mechanism, the 
GVM via replica symmetry breaking and the FRG via 
the generation of a cusp-like non-analyticity in R"(u), 
they are disconcertingly different in spirit. Physically, 
however, both capture the metastable states beyond the 
Larkin scale R c and should thus be related ||. A use- 
ful, quantitative and more general method, where this 
connection appears naturally, is still lacking. 

In this Letter we introduce such a method, encompass- 



ing both the FRG and the GVM. Via an exact calcula- 
tion of the effective action T[u] at large N, we obtain the 
FRG /3-function in any d at large N. Its detailed anal- 
ysis at dominant order, N — +oo, reveals that the FRG 
exactly reproduces (without invoking spontaneous RSB) 
the non-trivial result of the GVM for small overlap. The 
connections can be clarified using that T[u] also gives 
the probability distribution of a given mode u q . 0(1/N) 
corrections are computed, with the aim of understand- 
ing finite but large N. Further results, extensions and 
discussions will appear in p| . 

We start from the partition sum of an interface Zy = 
J T>[u] e~ n v[ u \/ T i n a given sample, with energy: 
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where J q = J w , 
mass m provides a scale. To obtain a non-trivial large N 
limit one defines the scaled field v = u/yN and chooses 
the distribution of the random potential O(N) rotation- 
ally invariant, e.g. its second cumulant as: 



V x (u)V x ,(u>)=R(u-u')6 xx , = NB((v-v') 2 )6 xx , (2) 



in terms of a function B(z) 
cumulants are scaled as 
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where the replica matrix field Xx = X x nas been in- 
troduced through a Lagrange multiplier. The bare in- 
teraction matrix potential U(x) — ^T?J2ab-B(Xab) — 
3mJ2abc S (Xab,Xbc,Xca) + ■■■ depends only on Xab = 
Xaa + Xbb — Xab — Xba and has & cumulant expansion in 
terms of sums with higher numbers of replicas. 

The effective action functional is defined as Legendre 
transform pi r [u] + W [J] = / J-u, with W[J] = In Z\j], 



J = V Nj. Its full calculation is given in [jllj. Since T[u] 
defines the renormalized 1PI vertices, its zero momen- 
tum limit defines the renormalized disorder. Thus we 
only need the result (per unit volume) for a uniform con- 
figuration of the replica field u" = u a — >/~N v a : 
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where vv stands for the matrix v a -Vb- We have computed 
the two first coefficients of the renormalized disorder in 
the 1/iV expansion U = U° + jjU 1 + . . .. Defining the 
notation d a bU {4>) = d^ ab U((f>) for any matrix </> with com- 
ponents 4> a b, we find at dominant order 



d ab u°(w) = d ab u(x(v)) 

X(v) =VV + T / [(q 2 
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i.e. a self consistent equation for dU°(vv), which, as we 
now show, contains both the GVM and the FRG. 

For simplicity, we now set all bare cumulants except 
B to zero. The above equations contain a huge amount 
of information, since they encode the full distribution 
(i.e. all cumulants) of the renormalized disorder, and 
are thus quite non-trivial to analyze. One limit where 
they "simplify" is when v is set to zero, since they then 
reproduce the Mezard Parisi (MP) equations || with 
x( w =0)a6 = J k G a b(k). These exhibit spontaneous RSB 
(with multiple solutions fl5f| ) and are solved by a hier- 
archical Parisi ansatz for x( v= ^)ab = x( w =0)(u) where 
< u < 1 is the overlap between replicas a and b. In 
the opposite limit of "strong" explicit symmetry breaking 
field (all v a b = v a — Vb ^ 0) we expect that the renormal- 
ized disorder U(vv) is given by a single saddle point and 
can be expanded in replica sums in terms of unambiguous 
renormalized cumulants, i.e. up to a constant: 
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This is the limit solved here, which we show is the one 
natural in the FRG, and amounts (in the RSB picture) to 
forcing the manifold in distant states. Work is in progress 
to analyze the rich crossover to RSB contained in @, 
when some of the v a b are set to zero. 

We can now expand, as detailed in fll|| , any quantity 
in (§) (e.g. a replica matrix M ab = M° b + £/ M\ hf + ... 
and its powers) in sums over an increasing number of 
free replica indices. This yields closed equations for the 
second cumulant (with I n := J k l/(fc 2 + m 2 )" ) 

B'{vlb) = B'{vlb + 2Th + Ah{B'{vl h ) - B'(0))) (7) 

with no other contributions from higher cumulants at any 
T. This is illustrated graphically in figure [l]. The three 
replica term S satisfies a closed equation involving only 
B, and all cumulants can be determined iteratively. 

The self consistent equation (Q) for the renormalized 
disorder can either be inverted directly (done below) or, 




FIG. 1: Top: typical T = contribution to B(v ab ). Bot- 
tom: self-consistent equation at leading order for B'(v^ b ) = 
B'(xab)- The wiggly line denotes a derivative, and is combi- 
natorially equivalent to choosing one B. At finite T one can 
attach an additional arbitrary number of tadpoles to any B. 



equivalently, turned into a FRG equation. We start with 
the solution B{x) — B(x) for m = oo (in presence of 
a fixed ultraviolet (UV) cutoff A) and then decrease m. 
Taking the derivative md m and rearranging gives: 

md m B'{x) = B"(x) [A(md m I 2 )(B'(x) - B'(0)) 

+2(md m TI 1 )(l + ^"(O))- 1 ] (8) 

valid for any d [JfTj. Since —^md rn Ii — m 2 l2, (||) has a 
well defined A — * +oo limit for d < 4. Then I 2 = A^ 22 ^- 
with A d = 2(47r)- d / 2 r[3 - d/2] and we can define the 
dimensionless function b(x) = 4Adrn A ''~ e B(xm~ 2 '>), ( for 
now arbitrary, which satisfies: 



-md m b(x) = (e- 4()b(x) + 2C,xb'(x) 

+\b'{x) 2 - b'(x)b'(0) + T m b'(x)(l 



(9) 

b"{Q)/e)- i 



where T m = T^m e , 6 = d-2 + 2(. The FRG equation 
that we have derived is valid, to dominant order in 1/N, 
in any dimension d < 4 and at any temperature T . Re- 
stricted to T = it correctly matches the one obtained 
by Balents and Fisher Q at any N but to lowest order in 
e = 4 — d. Furthermore, due to the self consistent equa- 
tion (0), it is fully integrable (not noted in Q). Indeed, 
(0) can be inverted into 
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where y = -b'(x), y = -b'(0) = -AA d m 2 ^ e B'(0), $ 
is the inverse function of —B'(x), i.e. (— B')(&(y)) = y 
and f m = 2Thm 2 ( (= T m /(2 - d) for d < 2). That this 
is also the general solution of the FRG equation can be 
seen by noting that (^) can be transformed into a linear 
equation for the inverse function x(y) 



md m x = (e - 2Qyx' + 2(x - (y - y ) + 
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with x' = x'(yo). The general T = solution of the ho- 
mogeneous part reproduces the $ term while a particular 
solution is x = {y — yo)/e using that md m yo = (2( — e)yo- 
To analyze the solutions of the large N FRG equation 
(^) , two approaches are legitimate, corresponding to dif- 
ferent points of view. The first, natural in mean field, 
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is exact integration. One discovers that ((^,(]9|) admit 
an analytic function as a solution, given by (|10|), only 
for to > m c where m r is the Larkin mass (more gen- 
erally a Larkin scale ||17)). Indeed its second derivative, 
^"(O)- 1 = B"(2TI 1 )- 1 -4I 2 , diverges (always for T = 0, 
d < 4, in some cases for T > 0) when to is lowered down 
to m c (which defines m c ). Since this expression is propor- 
tional to the replicon eigenvalue of the replica symmetric 
(RS) solution in the GVM, which exists for m > m c , 
the generation of a cusp in the FRG exactly coincides at 
large N with the instability of the RS solution. 

The second approach, natural in the RG Q , is to view 
the r.h.s. of (^|) as the large- iV limit of the true /3-function 
and to search for a zero. The general solution for 9 > 0, 
obtained from ([ll]), is parametrized by £: 

x = 7-2C + ^cT % y c (12) 

for ( > and ex = y - y - y \n(y/y ) for C = 0, with 
y = —b'(x). Here ya = —b'(0) is a fixed number. The 
value of the roughness exponent £ Q is selected by the 
decay of R(u) in (||) at large u, argued to be identical 
for B and B, i.e. if B'(z) ~ z -7 one finds C = C(t) = 
(4 — d)/2(l + 7) or £ = for shorter range correlations, 
to this order in 1/N. All the fixed points ( |l2|) have a 
cusp x'(yo) = (for £ ^ e/2) and are expected to be the 
physically correct solutions at small to. 

To show how to reconcile these two results, we study 
specific models, the long range (LR) correlations |J 
B(z) = i{l \ )Ad (a 2 + z) 1-7 and the short range (SR) 
gaussian correlator B(z) — jj-e~ z . Choosing £ = £(7) 
©yields: 

x = {y/~g)- Xh + e-\y - y ) - m 2C a 2 - f m (13) 

(and ln(gm~ e /y) in the r.h.s. for SR gaussian with_C = 0, 
a = 0). As one sees from Fig. ||, the r.h.s. of ( ^3[ ) has 
a minimum and decreasing m the curve x(y) cuts the 
axis x = closer to the minimum. It reaches it at 
m = m c where the solution acquires a cusp b'(0) — V(x) « 
y/-2(e-2()V(0)x and m 2 c <a 2 + f mc = {gi/e) 1/(1+l) = 
f c , 6'(0) = -g 1 /(i+7)( e / 7 )7/(i+7) ( = _ e f or SR with 
m e c = g/e). Although it is a priori unclear how to fol- 
low this solution for m < m c , the following remarkable 
property indicates how one may proceed. If we com- 
pute the /^-function, i.e. the r.h.s. of (g), using ( p^ ) at 
to = to c and C = £(7) we find that it exactly vanishes 
for all x > 0. Thus, for the potentials studied here, b(x) 
evolves according to ( |l3| ) until to c where it reaches its 
fixed point, and does not evolve for to < m c . This pro- 
vides unambiguously a solution beyond the Larkin scale, 
reconciles the two approaches and justifies the value ob- 
tained for (. The quantity B'^v 2 )- 1 = B"(x(w)) _1 -4I 2 
plays the role of a replicon eigenvalue and remains frozen 
and positive for v > 0. 

Here, for N = 00, temperature plays only a minor role 
in the case where disorder is relevant, i.e. for 0(7) > 




FIG. 2: The function x(y) given by (Q. 



(i.e. 2 < d < 4; d < 2 for 7 < 7 C = 2/(2 - d)), which 
is described by the T = fixed point ([12]). Contrarily 
to the one loop FRG, where b(x) remains analytic in a 
boundary layer x ~ T m at T > 0, here the denominator 
in (^) (resumming all orders in e) blows up and a cusp 
arises. Only in the marginal case (7 = 7 C ; d — 2 for 
SR disorder) we find a line of fixed points of (||) with 
Q = (2 — d)/2. For T > T c the disorder is analytic, given 
by (||) for all to (f m /f c = T/T c does not flow, and 
T c = 2tt in d = 2 for SR disorder). Below T c one recovers 
a cusp and the T = fixed point. Finally, for 7 > 7 C 
(d < 2) no cusp is generated as to — > and disorder is 
irrelevant, as in the corresponding RS solution of MP. 

Since the FRG aims at obtaining the universal be- 
haviour at small to, k of the correlation function of the 
manifold {v^v l L k ) = m d ^ 2{ 'g(k/m), we can now compute 
its zero momentum limit g(0), which is universal in the 
LR case, and compare with the GVM result #rsb(0). 
Extending || in presence of a mass, one shows that 
md m a m (u) = and thus the MP self energy is er m (u) = 
C u 2/0-i if u > u m , er m (u) = cr m (0) if u < u m . Within 
the FRG, <7frg (0) = -b'{0 + )/(2A d ) and substituting, we 
discover that 3frg(0) = T<r m (0)/m 4 < (?rsb(0), i.e. the 
FRG gives only but exactly the contribution from the 
overlaps u < u m (the distant states). This non-trivial 
information which within the MP approach requires a 
full RSB calculation, is obtained here without any RSB. 
Moreover using md m <j m (u) = one shows that contri- 
butions of larger overlaps in the MP result can be ob- 
tained by integrating the to dependent FRG result as 
5rsb(0) = T[a m (0)/m 4 + £° da m/ (0)/m' 4 +m- 2 - m - 2 }. 

These results can be understood as follows. The ef- 
fective action also describes the probability distribution 
Pv(w) in a given environment V, of the center of mass 
of the interface w — Nl } 2L d J x u x , i.e. one has r[{w a }] = 
— lirn^oo jfYp lnPy(wi) . . . Pv(w n ). Extension of the 
FRG beyond the Larkin scale requires giving a meaning 
to the u — + limit. We find here that what the FRG 
actually computes (from 6'(0 + )) is a second moment of 
w in presence of a small extra field V~Nj a such that all 
v a b 0; i- e - an average such that when there are sev- 
eral states the different replicas are chosen in maximally 
separated states (u = 0). 

In a previous study aiming to connect the RSB solution 
to the FRG |1 lnPy(u)) lnPy(— w) was computed and 
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FIG. 3: Contribution to the second cumulant at order l/N. 



used to define a renormalized second cumulant of the 
disorder. This quantity is however different from the one 
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(3(b) = eb+^b' 2 -b'b'v 



in the FRG, obtained here, and does not reproduce the 
second moment of w, neither <?(0)r.sb, nor 3(0)frg- In 
addition, since ||] used the unperturbed MP saddle point, 
the two calculations focus on different regimes (i£ 
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in progress to connect these regimes, and obtain a more 
complete version of the FRG, using our equations (||) and 
summing over RSB saddle points [[15| . 

The FRG approach should allow a quantitative study 
of finite N beyond possible artifacts of N — oo, and 
extension of Q to any d. The calculation of the correction 
at order 1/iV to B is involved, and was performed using 
two complementary methods, a graphical one, see Fig. |^, 
and the algebraic expansion in number of replica sums. 
The resulting expression for the /3-function at order l/N 
and T = is UV-convergent and reads for A — > oo 

2exlf(b' -b>)b"(l-I 2 b")H- 2 



+eI 4 p (&; ) -60V(2+(2/ 2 -/f)fo'0ff- 2 -^^ 

H. 



-2lE 2 (b' -brb''H- 1 +2ex 2 I^(l+(I 2 +ll)b'')H- 



l+(/ 2 -/f)6", 



V 



(14) 



with I 2 = e \ Jf = Jf 1; if = Jf 2 , if = Jf 2 and J^ m = 

-hI k (^) n (Io +P r, b< = b'(x), b' = b'(0), etc... The 
corresponding expressions at T > have been obtained. 
Analysis of these formidable expressions is in progress. 
In particular, we have not included in (|l4|) anomalous 
terms arising from the non-analytic structure. We have 
checked that (|l4| ) non-trivially reproduces the two loop 
FRG equation for the A^-component model of JlO|, pT[ . 

In summary, via an exact calculation of the effective 
action at large N we have derived equations valid in any 
d containing both the GVM and the FRG. The FRG 
and its continuation to m < m c are consistent with the 
main results of the full and the marginal one step RSB- 
solutions of MP. Since it reproduces the non-trivial small 



overlap results it provides another way to attack finite N. 
l/N corrections have been obtained. Our study hints at 
further connections between: l/N and thermal bound- 
ary layers, RSB, and how to fix the ambiguities in the 
anomalous terms in the /^-function. Their understanding 
should allow quantitative progress in the SR case (e.g. for 
d — 1, equivalent to the A^-dimensional KPZ equation). 
Applications of the method to other complex systems is 
in progress. We have also computed the effective action 
for the random field 0(Af)-model at large A^ p] . Fi- 
nally, applications to the dynamics offer the hope to go 
systematically beyond mode-coupling approximations. 
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